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ABSTRACT 



We report on the 7-ray activity of the high-synchrotron-peaked BL Lacer- 
tae object Mrk421 during the first 1.5 years of Fermi operation, from 2008 
August 5 to 2010 March 12. We find that the Large Area Telescope (LAT) 
7-ray spectrum above 0.3 GeV can be well-described by a power-law function 
with photon index T = 1.78 ± 0.02 and average photon flux F(> 0.3 GeV) = 
(7.23 ± 0.16) x 10 -8 phcm~ 2 s _1 . Over this time period, the Fermi-LAT spec- 
trum above 0.3 GeV was evaluated on 7-day-long time intervals, showing sig- 
nificant variations in the photon flux (up to a factor ~ 3 from the minimum 
to the maximum flux), but mild spectral variations. The variability amplitude 
at X-ray frequencies measured by RXTE/ASM and Swift/BAT is substantially 
larger than that in 7-rays measured by Fermi-LAT, and these two energy ranges 
are not significantly correlated. We also present the first results from the 4.5- 
month-long multifrequency campaign on Mrk421, which included the VLBA, 
Swift, RXTE, MAGIC, the F-GAMMA, GASP-WEBT, and other collaborations 
and instruments which provided excellent temporal and energy coverage of the 
source throughout the entire campaign (2009 January 19 to 2009 June 1). During 
this campaign, Mrk421 showed a low activity at all wavebands. The extensive 
multi-instrument (radio to TeV) data set provides an unprecedented, complete 
look at the quiescent spectral energy distribution (SED) for this source. The 
broad band SED was reproduced with a leptonic (one-zone Synchrotron Self- 
Compton) and a hadronic model (Synchrotron Proton Blazar). Both frameworks 
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are able to describe the average SED reasonably well, implying comparable jet 
powers but very different characteristics for the blazar emission site. 

Subject headings: acceleration of particles — radiation mechanisms: non-thermal 
- galaxies: active — BL Lacertae objects: general — BL Lacertae objects: 
individual (Mrk421) — gamma rays: observations 



1. Introduction 



Blazars are active galaxies believed to have pairs of relativistic jets flowing in opposite 
directions closely aligned to our line of sight. Their spectral energy distributions (SEDs) are 
dominated by beamed jet emission and take the form of two broad nonthermal components, 
one at low energy, peaking in the radio through optical, and one at high energies, peaking 
in the 7-rays. Some blazars have been well-monitored for decades and along a wide range of 
wavelengths. Although there is ample evidence for the electron synchrotron origin of the low- 
energy bump, the existing data do not allow an unambiguous identification of the radiation 
mechanism responsible for the high-energy bump. One reason for this is that the high-energy 
bump is poorly constrained due to the lack of observations at energies between ~ 0.1 MeV 
and 0.3 Te V. This gap was filled to some extent by EGRET on the Compton Gamma-Ray 
Observatory (jHartman et al.lll999l ). However, its moderate sensitivity and limited observing 
time precluded detailed cross-correlation studies between 7-ray and lower-energy wavebands. 
On the other hand, the current generation of TeV imaging atmospheric Cherenkov telescopes 
(IACTs), HESS, MAGIC, and VERITAS, which have good sensitivity at energies as low as 
0.1 TeV, did not start scientific operation until 2004; that is, well after EGRET had stopped 
operating. 

This has changed with the launch of the Fermi Gamma-ra y Space Telescope in June 
2008. In science operation since 2008 August, its LAT instrument (lAtwood et al.ll2009l ) views 
the entire sky in the 20 MeV to greater than 300 G eV range every three hours. The one 
year LAT AGN Catalog (1LAC; lAbdo et al.l l2010bl ) contains around 600 blazars, a factor 
of ~ 10 greater than EGRET detected during its entire operational lifetime. For the first 
time, simultaneous observations of Fermi with the latest generation of IACTs can cover the 
entire high-energy bump. Combining this with simultaneous low-energy observations gives 
an unprecedented multiwavelength view of these enigmatic objects. 

Blazars found in low states are particularly poorly studied. This is due in part to the 
lower sensitivity of previous instruments, and in part to the fact that multiwavelength mon- 
itoring programs, including space-based instruments, are mostly triggered when an object 
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enters a particularly bright state, as observed by gr ound-based optical telescopes and all- 
sky monitors such a s the RXTE faradt et al.l Il993h All Sky Monitor (ASM) or the Swift 
( IGehrels et al.l 12004 ) Burst Alert Telescope (BAT). Having a well-measured low-state SED 
will be useful for constraining models and as a baseline to which other, flaring states can be 
compared. This will be crucial for answering many of the questions regarding these objects. 

Markarian 421 (Mrk421; RA=ll h 4 m 27.31 s , Dec=38° 12' 31.8" , J2000, redshift 
z = 0.031) is a high-synchrotr on-peaked (HSP) BL Lac (according to the classification 
presented in lAbdo et al.l ( 12010d )). that is one of the brightest sources in the extragalactic 
X-ray/TeV s ky. Mrk421 was a ctually the first extragalactic object to be di scovered as a 



TeV emitter (IPunch et al.lll992l ). and one of the fastest varying 7-ray sources (IGaidos et al. 



19961 ). During the last decade, there were a large number of publications on the very high 
energy (VHE) 7-ray spectrum of this source, which has been measured with almost all the ex- 
isting IACTs fcreimrich et al.l l2002[ kharonian et al.ll2002L I2OO3L I2OO5I : I Albert et al1l2007al : 
Acciari et al.l 120091 ) . Among other things, we learned that the source shows evidence for a 
spectral hardening with increasing flux. The SED and the multifrequency correlations of 
Mrk421 have also been int ensively studied in the past through dedicated multifrequency 



obseryations of the sourc e (IKatarzvhs 



2006: Fossati et al. 2008: Horan et al. 



ri et al.l l2003t iBlazejowski et al.l 120051 : iRevillot et al. 



2009), which showed a positive but very complex 



relation between X-rays and VHE 7-rays, and that a simple one-zone Synchrotron Self- 
Compton model (SSC) with an electron distribution parameterized with one or two power- 
laws seemed to describe the collected SED well during the observing campaigns. During a 
strong flare in June 2008, the source was also detected with the gamma-ray telescope AGILE 
an d, for the first time, a hint of correlation between optical and TeV energies was reported 
by iDonnarumma et al.l (120091 ) . 



Despite the large number of publications on Mrk421, the details of the physical processes 
underlying the blazar emission are still unknown. The main reasons for this are the sparse 
multifrequency data during long periods of time, and the moderate sensitivity available in the 
past to study the 7-ray emission of this source. In addition, as occurs often with studies of 
blazars, many of the previous multifrequency campaigns were triggered by an enhanced flux 
level at X-rays and/or 7-rays, and hence many of the previous studies of this source are biased 
towards "high-activity" states, where perhaps distinct physical processes play a dominant 
role. Moreover, we have very little information from the MeV-GeV energy range: 9 years of 
operation with EGRET resulted in only a f ew viewing periods w ith a signal significance of 
barely 5 standard deviations (a, hereafter) ( iHartman et al.lll999l ). which precluded detailed 
correlation studies with other energy bands. 



We took advantage of the new capabilities provided by Fermi-LAT and the new IACTs, 
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as well as the existing capabilities for observing at X-ray and lower frequencies, and orga- 
nized a multifrequency (from radio to TeV) campaign to observe Mrk 421 over 4.5 months. 
The observational goal for this campaign was to sample Mrk 421 every 2 days, which was 
accomplished at optical, X-rays and TeV energies whenever the weather and/or technical 
operations allowed. The Fermi-LAT operated in survey mode and thus the source was con- 
stantly observed at 7-ray energies. In this paper, we report the overall SED averaged over 
the duration of the observing campaign. A more in-depth analysis of the multifrequency 
data set (variability, correlations and implications) will be given in a forthcoming paper. 

The work is organized as follows: In section 2 we introduce the LAT instrument and 
report on the data analysis. In section 3 we report the flux/spectral variability in the 7- 
ray range observed by Fermi-LAT during the first 1.5 years of operation, and compare it 
with the flux variability obtained with RXTE/ASM and Swift /BAT, which are also all-sky 
instruments. In section 4 we report on the spectrum of Mrk 421 measured by Fermi, and 
section 5 reports on the overall SED collected during the the 4.5 month long mult i- wavelength 
campaign organized in 2009. Section 6 is devoted to SED modeling of the multifrequency 
data with both a hadronic and a leptonic model, and in section 7 we discuss the implications 
of the experimental and modeling results. Finally, we conclude on section 8. 



2. Fermi-LAT Data Selection and Analysis 

The Fermi-LAT is a 7-ray telescope operating from 20 MeV to > 300 GeV. The instru- 
ment is an array of 4 x 4 identical towers, each one consisting of a tracker (where the photons 
are pair-converted) and a calorimeter (where the energies of the pair-converted photons are 
measured). The entire instrument is covered with an anticoincidence detector to reject the 
charged-particle background. The LAT has a large peak effective area (0.8 m 2 for 1 GeV 
photons), an energy resolution typically better than 10%, and a field of view (FoV) of about 
2.4 sr with an angular resolution (68% containment angle) bet ter than 1° for energi es above 



1 GeV. Further details on the description of LAT are given by lAtwood et al.l (120091 ) . 



The LAT data reported in this paper were collected from 2008 August 5 (MJD 54683) to 
2010 March 12 (MJD 55248). During this time, the Fermi-LAT instrument operated almost 
entirely in survey mode. The analysis was performed with the Science Tools software package 
version v9rl5p6. Only events having the highest probability of being photons, those in the 
"diffuse" class, were used. The LAT data were extracted from a circular region with a 10° 
radius centered at the location of Mrk 421. The spectral fits were performed using photon 
energies greater than 0.3 GeV, where the effective area of the instrument is large (> 0.5 m 2 ) 
and the angular resolution relatively good (68% containment angle smaller than 2°). The 



- 11 - 



spectral fits using energies above 0.3 GeV are less sensitive to possible contamination from 
non-accounted (transient) neighboring sources, and have smaller systematic errors, at the 
expense of reducing somewhat the number of photons from the source. In addition, a cut 
on the zenith angle (< 105°) was also applied to reduce contamination from the Earth limb 
7-rays, which are produced by cosmic rays interacting with the upper atmosphere. 

The background model used to extract the 7-ray signal includes a Galactic diffuse emis- 
sion component and an isotropic component. The model that we adopted for the Galactic 
component is given by the file gll_iem_v02.fit, and the isotropic component, which is the 
sum of the extragalactic diffuse emission and the residual charged particle background, is 



parametrized by the file isotropic_iem_v02 l 115 l . The normalization of both components in 



the background model were allowed to vary freely during the spectral point fitting. The 
spectral analyses (from which we derived spectral fits and photon fluxes) were 
performed with the post-launch instrument response functions P6_V3_DIFFUSE us- 
ing an unbinned maximum likelihood method. The systematic uncertainty in 
the flux were estimated as 10% at 0.1 GeV, 5% at 560 MeV and 20% at 10 GeV and 
above 



116 



3. Flux and Spectral Variability 

The sensitivity of Fermi-LAT is sufficient to accurately monitor the 7-ray flux of Mrk 421 



on short timescales (few days)l 117 l. The measured 7-ray flux above 0.3 GeV and the photon 
index from a power law (PL) fit are shown in Figure __( The data span the time from 2008 
August 5 (MJD 54683) to 2009 March 12 (MJD 55248) and they are binned on time intervals 



of 7 days. The Test Statistic (TS) valued 118 ! for the 81 time intervals are typically well in 



excess of 100 (~10 sigma). The number of intervals with TS<100 is only 9 (11%). The 
lowest TS value is 30, which occurs for the time interval MJD 54899-54906. This low signal 
significance is due to the fact that the Fermi-LAT instrument did not operate during the 



http : //f ermi . gsfc.nasa.gov/ssc/data/access/lat/BackgroundModels .html 



115 

116 See http : //f ermi . gsf c . nasa.gov/ssc/data/analysis/LAT_caveats .html 

117 The number of photons from Mrk 421 (above 0.3 GeV) detected by LAT in one day is 
typically about 6. 

118 The Test Statistic TS=2Alog(likelihood) between models with and without the source is a measure 
of the probability of havi ng a point 7-ray so urce at the location specified. The TS value is related to the 



significance of the signal ( Mattox et al. 19961) 
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time interval MJD 54901-5490a 119 l and hence only three out of the seven days of the interval 
contain data. The second lowest TS value is 40, which occurred for the time interval 54962- 
54969. During the first 19 months of Fermi operation, Mrk421 showed relatively mild 7-ray 
flux variations, with the lowest photon flux F(> 0.3 GeV) = (2.6 ± 0.9) x 10~ 8 cm~ 2 s _1 
(MJD 54906-54913; TS=53) and the highest F{> 0.3 GeV) = (13.2 ± 1.9) x lO^cm-V 1 
(MJD 55200-55207; TS=355). A constant fit to the flux points from Figure GDgave a \ 2 = 159 
for 82 degrees of freedom (probability the flux was constant is 8 x 10~ 7 ), hence indicating 
the existence of statistically significant flux variability. On the other hand, the photon index 
measured in 7-day- long time intervals is statistically compatible with being constant, as 
indicated by the results of the constant fit to all the photon index values, which gave x 2= 
87 for 82 degrees of freedom (probability no variability is 0.34). The scatter plot with Flux 
versus Index in Figure [H shows that there is n o obvious relation between these two quantities. 



We quantified the correlation as prescribed in lEdelson &: Krolikl (119881 ). obtaining a Discrete 



Correlation Function DCF = 0.06 ± 0.11 for a time lag of zero. 

It is interesting to compare the 7-ray fluxes measured by Fermi with those historical 



ones recorded by EGRET. From the 3rd EGRET catalog (IHartman et al.lll999l ). one can 
see that the highest /lowest significantly- measured (TS > 25) photon fluxes are F Max (> 0.1 
GeV) = (27.1±6.9)xl0~ 8 cm- 2 s" 1 (TS=32) and F Min (> 0.1 GeV) = (10.9±2.8) x HT 8 cm" 2 
s- 1 (TS=26); where F(> 0.1 GeV) is the flux above 0.1 GeV. These values do not deviate 
by more than 2 sigmas from the P1234 average, F(> 0.1 GeV) = (13.8 ± 1.8) x 10 8 
cm -2 s 1 (TS=100), and hence EGRET did not detect significant variability in 
the flux from Mrk421. We can easily obtain the Fermi F(> 0.1 GeV) fluxes by using 
the flux (F) index (T) values reported in Figure Q] (E>0.3 GeV): F{> 0.1 GeV) = F(> 0.3 
GeV) x (0.3/0. l) ri . Applying this simple formalism one gets, for the max/min fluxes from 
Figure [fl F Max (> 0.1 GeV) = (25.7 ± 4.7) x 10~ 8 cm~ 2 s" 1 and F Mm (> 0.1 GeV) = 
(5.6 ±2.4) x 10~ 8 cm 2 s~ x . The maximum flux measured by EGRET and LAT are 
similar, although the minimum fluxes are not. The LAT's larger effective area 
compared to EGRET permits detection of lower 7-ray fluxes. In any case, the 
EGRET and LAT fluxes are comparable, which may indicate that Mrk421 is not 
as variable in the M e V / GeV range as at other wavelengths, particularly X-rays 



and TeV 7-rays (e.g. IWagnerl 12008) . 



The Fermi-LAT capability for constant source monitoring is nicely complemented at 
X-ray energies by RXTE/ASM and Swift /BAT, the two other all-sky instruments which 
can probe the X-ray activity of Mrk421 on 7-day-long time intervals. Figure [2] shows the 



119 



The LAT did not operate during the time interval MJD 54901-54905 due to an unscheduled shutdown. 
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measured fluxes by ASM in the energy range 2 — lOkeV, by BAT at 15 — 50keV, an d by 
LAT in two different energy bands: 0.2 — 2 GeV (low energy) and > 2 GeV (high energy' 



120 



The low/high Fermi-LAT energy bands were chosen (among other reasons) to 
produce comparable flux errors. This might seem surprising at first glance, 
given that the number of detected photons in the low energy band is about 5 
times larger than in the high energy band (for a differential energy spectrum 
parameterized by a power law with photon index of 1.8, which is the case of 
Mrk421). Hence the number of detected 7-rays decreases from about 50 down 
to about 10 for time intervals of 7 days. The main reason for having comparable 
flux errors in these two energy bands is that the diffuse background, which follows 
a power law with index 2.4 for the high galactic latitude of Mrk421, is about 25 
times smaller in the high energy band. Consequently, Signal/Noise ~ Ng/ ^/(Nb) 
remains approximately equal. 

We do not see variations in the LAT hardness ratio (i.e. F(> 2GeV)/F(0.2 — 
2 GeV) with the 7-ray flux, but this is limited by the relatively large uncertain- 
ties and the low 7-ray flux variability durin g th is time interval. The data from 



RXTE/ASM were obtained from the ASM web pagd 121 !. We filtered out the data according 
to the provided prescription in the ASM web page, and made a weighted average of all the 
dwells (scan/rotation of the ASM Scanning Shadow Cameras lasting 90 seconds) from the 
7-day- long time intervals defined for the Fermi data. The data from Swift /BAT were gath- 



ered from the BAT web page I. We retrieved the daily averaged BAT values and produced 



weighted average for all the 7-day-long time intervals defined for the Fermi data. 

The X-ray flux from Mrk 421 was ~ 1.7 ct s 1 in ASM and ~ 1.9 x 10" 3 
ct s -1 cm -2 in BAT. These fluxes correspond to ~22 mCrab in ASM (1 Crab = 
75 ct s x ) and 9 mCrab in BAT (1 Crab = 0.22 ct s _1 cm~ 2 ), although given the 



recent reports on flux va r iability from the C rab Nebula (see Wilson-Hodge et al. 



201 ll : lAbdo et al.ll2011al : iTavani et al.ll201lh . the flux from the Crab Nebula is 



not a good absolute standard candle any longer and hence those numbers need 
to be taken with caveats. One may note that the X-ray activity was rather low during 
the first year of Fermi operation. The X-ray activity increased around MJD 54990 and then 



120 The fluxes depicted in the light curves were computed fixing the photon index to 1.78 (average index 
during the first 1.5 years of Fermi operation) and fitting only the normalization factor of the power law 
function. 

121 See http : //xte .mit . edu/ASM_lc .html 

122 See http : / / swift . gsf c . nasa.gov/docs/ swift/results/transients/ 
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increased even more around MJD 55110. The 7-ray activity seemed to follow some of the 
X-ray activity, but the variations in the 7-ray range are substantially smoother than those 
observed in X-rays. 

Figure [3] shows the same light curves as Figure [21 but only during the period of time 
after MJD 55110 (when Mrk421 showed high X-ray activity) with a time bin of only 3 
days. During this time period the ASM and BAT flux (integrated over 3 days) went beyond 
5 ct s _1 and 8 x 10~ 3 ct s _1 cm -2 , respectively, which implies a flux increase by a factor of 
5-8 with respect to the average fluxes during the first year. It is worth noting that these 
large flux variations do not have a counterpart at 7-ray energies measured by Fermi-LAT . 
The MeV/GeV flux measured by LAT remained roughly constant, with the exception of a 
flux increase by a factor of ~2 for the time interval around MJD 55180-55210 and around 
MJD 55240-55250, which was also seen by RXTE/ASM and (to some extent) by Swift /BAT. 

We quantified the correlatio n among the light curves shown in Figure [2] and Figure 



[3] following the prescription from lEdelson fc Krolikl (119881 ). The results are shown in 



table [JJ for a time lag of zero, which is the one giving the largest DCF values. 
There is no indication of correlated activity at positive/negative time lags in 
the DCF vs time plot for any of the used X-ray /7-ray bands. The advantage of 
using the DCF instead of the Pearson's correlation coefficient is that the later 
does not consider the error in the individual flux points, while the former does. 
In this particular situation it is relevant to consider these errors because they 
are sometimes comparable to the magnitude of the measured flux variations. 
The main result is a clear (DCF / DCF error ~ 4) correlation between ASM and 
BAT, while there is no indication of X-ray /7-ray correlation (DCF/DCF error < 2). 
The correlation between the Fermi-LAT fluxes below and above 2 GeV is not 
significant (DCF = 0.31 ± 0.14), which is probably due to the low variability at 
7-rays, together with the relatively large flux errors for the individual 7-day-long 
and 3-day-long time intervals. 



We followed the prescription given in IVaughan et all (120031 ) to quantify the flux vari- 



ability by means of the fractional variability parameter F var , as a function of energy. In order 
to account for the individual flux measurement errors (cr prr A we used the "excess variance" 



(INandra et all 119971 ; lEdelson et all |2002| ) as an estimator of the intrinsic source variance. 
This is the variance after subtracting the expected contribution from measurement errors. 
For a given energy range, the F var is calculated as 



S 2 <c ^"g rr ^ 



< F > 2 

where < F > is the mean photon flux, S the standard deviation of the N flux points, and 
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< <7g rr > the average mean square error, all determined for a given energy bin. 

Figure H] shows the derived F var values for the 4 energy ranges and time window covered 
by the light curves shown in Figure [2J The fractional variability is significant for all energy 
ranges, with the X-rays having a substantially higher variability than the 7-rays. 

It is interesting to note that, while the PL photon index variations from Figure[T]were not 
statistically significant (x 2 /NDF = 87/82), Figure H] shows that the fractional variability 
for photon energies above 2 GeV is higher than that below 2 GeV; specifically F var (E < 
2 GeV) = 0.16 ± 0.04 while F var (E > 2GeV) = 0.33 ± 0.04. This apparent discrepancy 
between the results reported in Figure [1] and the ones reported in Figure 0] (produced with 
the flux points from Figure [2]) might be due to the fact that, on time scales of 7-days, the 
photons below 2 GeV dominate the determination of the PL photon index in the unbinned 
likelihood fit. In other words, the source is bright enough in the energy range 0.3 — 2 GeV 
such that the (relatively few) photons above 2 GeV do not have a large (statistical) weight 
in the computation of the PL photon index. Consequently, we are more sensitive to spectral 
variations when doing the analysis separately for these two energy ranges. 

One may also note that, beside the larger variability in the Fermi fluxes above 2 GeV 
with respect to those below 2 GeV, the variability in the BAT fluxes (15-50 keV) is also 
higher than that of ASM (2-10 keV). The implications of this experimental result will be 
further discussed in section 17.31 in light of the modeling results presented in section 16.21 

4. Spectral Analysis up to 400 GeV 

The LAT instrument allows one to accurately reconstruct the photon energy over many 
orders of magnitude. Figure [5] shows the spectrum of Mrk421 in the energy range 0.1- 
400 GeV. This is the first time that the spectrum of Mrk421 can be studied with this 
level of detail over this large fraction of the electromagnetic spectrum, which includes the 
previously unexplored energy range 10-100 GeV. The spectrum was computed using the 
analysis procedures described in section [2j In order to reduced systematics, the spectral fit 
was performed using photon energies greater than 0.3 GeV, where the LAT instrument has 
good angular resolution and large effective area. The black line in Figure [5] is the result 
of a fit with a single PL function over the energy range 0.3-400 GeV, and the red contour 
is the 68% uncertainty of the fit. The data are consistent with a pure PL function with 
photon index of 1.78 ± 0.02. The black data points are the result of performing the analysis 
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on differential energy ranges (2.5 bins per decade of energ y)H The points are well within 



l-2er from the fit to the overall spectrum (black line), which confirms that the entire Fermi 
spectrum is consistent with a pure power law function. 

However, it is worth noticing that the error bars at the highest energies are relatively 
large due to the low photon count. In the energy bins 60-160 GeV and 160-400 GeV, the 
predicted (by the model for Mrk421) number of photons detected by LAT is 33 and 11, 
respectively. Even though the low background makes those signals very significant (TS 
values of 562 and 195 respectively), the statistical uncertainties in the energy flux values are 
naturally large and hence they could hide a potential turnover in the spectrum of Mrk421 
at around 100 GeV. Indeed, when performing the likelihood analysis on LAT data above 
100 GeV, one obtains a photon flux above 100 GeV of (5.6 ±1.1) x 10 -10 phcm -2 s _1 with 
a photon index of 2.6 ± 0.6, which might suggest a turnover in the spectrum, consistent 



with the TeV spectra determined by past observations with IACTs (IKrennrich et al.l 12002 



Aharonian et al.ll2003l . l2005t Albert et al.ll2007ai ). In order to make a statistical evaluation 
of this possibility, the LAT spectrum (in the range 0.3-400 GeV) was fit with a broken power 
law (BPL) function, obtaining the indices of 1.77 ± 0.02 and 2.9 ± 1.0 below and above the 
break energy of 182 ±39 GeV, respectively. The likelihood ratio of the BPL and the PL gave 
0.7, which, given the 2 extra degrees of freedom for the BPL function, indicates that the 
BPL function is not statistically preferred over the PL function. Therefore, the statistical 
significance of the LAT data above 100 GeV is not sufficiently high to evaluate the potential 
existence of a break (peak) in the spectrum. 



5. Spectral Energy Distribution of Mrk421 during the 4.5 month long 

Multifrequency Campaign from 2009 

As mentioned in ^TJ we organized a multifrequency (from radio to TeV photon energies) 
campaign to monitor Mrk421 during a time period of 4.5 months. The observing campaign 
started on 2009 January 19 (MJD 54850) and finished on 2009 June 1 (MJD 54983). The 
observing strategy for this campaign was to sample the broad-band emission of Mrk421 
every 2 days, which was accomplished at optical, X-ray and TeV energies when the weather 
and/or technical limitations allowed. The main goal of this project was to collect an extensive 
multifrequency data set that is simultaneous and representative of the average/typical SED 



123 Because of the analysis being carried out in small energy ranges, we fixed the spectral index to 1.78, 
which is the value obtained when fitting the entire energy range. We repeated the same procedure fixing the 
photon indices to 1.5 and 2.0, and found no significant change. Therefore, the results from the differential 
energy analysis are not sensitive to the selected photon index used in the analysis. 
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from Mrk421. Such a data set can provide additional constraints that will allow us to refine 
the emission models, which in turn will provide new insights into the processes related to 
the particle acceleration and radiation in this source. In this section we describe the source 
coverage during the campaign, the data analysis for several of the participating instruments, 
and finally we report on the averaged SED resulting from the whole campaign. 



5.1. Details of the Campaign: Participating Instruments and Temporal 

Coverage 

The list of all the instruments that participated in the campaign are reported in Ta- 



ble|2j and the scheduled observations can be found online 124 !. We note that in some cases 
the planned observations could not be performed due to bad observing conditions, while on 
other occasions the observations were performed but the data could not be properly analyzed 
due to technical problems or rapidly changing weather conditions. Figure^ shows the time 
coverage as a function of the energy range for the instruments/observations used to pro- 
duce the SED shown in Figure[8l Apart from the unprecedented energy coverage (including, 
for the first time, the GeV energy range from Fermi-LAT), the source was sampled very 
uniformly with the various instruments participating in the campaign and, consequently, it 
is reasonable to consider the SED constructed below as the actual average (typical) SED 
of Mrk421 during the time interval covered by this multifrequency campaign. The largest 
non-uniformity in the sampling of the source comes from the Cherenkov Telescopes, which 
are the instruments most sensitive to weather conditions. Moreover, while there are many 
radio/optical instruments spread all over the globe, in this observing campaign, only two 
Cherenkov Telescope observatories participated, namely MAGIC and Whipple. Hence, the 
impact of observing conditions was more important to the coverage at the VHE 7-ray ener- 
gies. During the time interval MJD54901-54905 the Fermi satellite did not operate due to a 
spacecraft technical problem. The lack of Fermi-LAT data during this period is clearly seen 
in FigureEl 

We note that Figure[6] shows the MAGIC and Whipple coverage in VHE 7-ray energies, 
but only the MAGIC observations were used to produce the spectra shown in FigureEJ The 
more extensive, but less sensitive, Whipple data (shown as gr ey boxes in FigureE]) were 



primarily taken to determine the light curve (jPichel et all 120091 ) and a re-optimization was 



required to derive the spectrum, which will be reported elsewhere. 



124 See https : //confluence . slac . Stanford . edu/display/GLAMCDG/Campaign+on+Mrk421+ ( Jan+2009+to+May+2009) ) 

maintained by D. Paneque. 
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In the following paragraphs we briefly discuss the procedures used in the analysis of 
the instruments participating in the campaign. The analysis of the Fermi-LAT data was 
described in §2] and the results obtained will be described in detail in §5.21 



5.1.1. Radio Instruments 



Radio data were taken for this campaign from single-dish telescopes, one mm-interferometer, 
and one Very Long Baseline Interferometry (VLBI) array, at frequencies between 2.6 GHz 
and 225 GHz (see Table[5]). The single-dish telescopes were the Effelsberg 100 m radio tele- 
scope, the 32 m Medicina radio telescope, the 14 m Metsahovi radio telescope, the 32 m Noto 
radio telescope, the Owens Valley Radio Observatory (OVRO) 40 m telescope, the 26 m 
University of Michigan Radio Astronomy Observatory (UMRAO) and the 600 meter ring 
radio telescope RATAN-600. The mm-interferometer was the Sub-millimeter Array (SMA). 
The NRAO Very Long Baseline Array (VLBA) was used for the VLBI observations. For 
the single-dish instruments and SMA, Mrk421 is pointlike and unresolved at all observing 
frequencies. Consequently, the single-dish measurements denote the total flux density of 
the source integrated over the whole source extension . Details of the observ ing strategy 



and data reduction are given bvlFuhrmann et & 



projec t), iTerasranta et al 



19981. Metsah o vi), 



(200 81): lAngelakis et all (120081. F-GAMMA 



Aller et all (119851 . UMRAO), 



Kovalev et all ( 119991 . RATAN-600) and lRichards et all (12011 in 



Venturi et al. 



(120011 . Medicina and Noto). 
preparation, OVRO). 

The VLBA data were obtained at various frequencies (5, 8, 15, 24, 43 GHz) through 
various programs (BP143, BK150 and MOJAVE). The data were redu ced following stan 



dard procedures for d ata reduction and calibration (see, for example, iLister et all 12009 



Sokolovsky et al.ll2010l . for a description of the MOJAVE and BK150 programs, respectively). 
Since the VLBA angular resolution is smaller than the radio source extension, measurements 
were performed for the most compact core region, as well as for the total radio structure 
at parsec scales. The core is partially-resolved by our 15, 24 and 43 GHz observat i ons ac - 
cording to the resolution criterion proposed by iKovalev et al.1 (120051 ) and lLobanovl (120051 ). 
The VLBA core size was determined with two-dimensional Gaussian fits to the measured 
visibilities. The FWHM size of the core was estimated to be in the range 0.06-0.12 mas at 
the highest observing frequencies, 15, 24 and 43 GHz. Both the total and the core radio flux 
densities from the VLBA data are shown in FigureEJ 
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5.1.2. Optical and Near-IR Instruments 



The coverage at optical frequencies was provided by various telescopes around the globe, 
and this decreased the sensitivity to weather and technical difficulties and provided good 
overall coverage of the source, as depicted in F igureEl Many of the obs ervations were per- 



formed within the GASP-WEBT program (e.g.. IVillata et al.ll2008l . 120091 ); this is the case for 
the data collected by the telescopes at Abastumani, Lulin, Roque de los Muchachos (KVA), 
St. Petersburg, Talmassons, and Valle d'Aosta observatories (R band). In addition, the tele- 
scopes GRT, ROVOR, New Mexico Skies and MitSume provided data with various optical 
filters, while OAGH and WIRO provided data at near- infrared wavelengths. See Tabled for 
further details. 



All the optical and near/IR instruments used the calibration stars reported inlVillata et al. 



( 1998 ). and the Galactic extinction was corrected with the coefficients given in lSchlegel et al. 
( 119981 ). The flux from the host galaxy (which is significant only below v ~ 10 15 Hz) was esti- 



m ated using the flux val ues at the R band from iNilsson et al.l ( 120071 ) and the colors reported 



m 



Fukugita et al.l ( 119951 ) . and then subtracted from the measured flux. 



5.1.3. Swift/UVOT 



The Swift Ultraviolet /Optical Telescope (UVOT; iRoming et al.l 120051 ) dataset includes 
all the observations performed during the time interval MJD 54858 to 54979, which amounts 
to 46 single pointing observations that were requested to provide UV coverage during the 
Mrk421 multifrequency campaign. The UVOT telescope cycled through each of three ul- 
traviolet passbands (UVW1, UVM2, UVW2). Photometry was computed using a 5arcsec 
source region around Mrk 421 us ing a custom UVOT pipeline that performs the calibrations 
presented in iPoole et al.l (120081 ). Moreover, the custom pipeline also allows for separate 



observation-by-observation, corrections for astrometric mis-alignments (lAcciari et al.l 12010 
in preparation). A visual inspection was also performed on each of the observations to ensure 
proper data quality selection and correction. The flux measurements o btained have been cor- 
rected for Galacti c extinction E B _ V = 0.019 mag ( ISchlegel et al.lll9981 ) in each spectral band 
dFitzpatricklll999l ). 



5.1.4. Swift/XRT 



All the Swift X-ray Telescope (XRT; iBurrows et al.l 120051 ) Windowed Timing observa- 
tions of Mrk 421 carried out from MJD 54858 to 54979 were used for the analysis: this 
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amounts to a total of 46 observations that were performed within this dedicated multi- 
instrument effort. The XRT data set was first processed with the XRTDAS software package 
(v.2.5.0) developed at the ASI Science Data Center (ASDC) and distributed by HEASARC 
within the HEASoft package (v. 6. 7). Event files were calibrated and cleaned with standard 
filtering criteria with the xrtpipeline task using the latest calibration files available in the 
Swift CALDB. The individual XRT event files were then merged together using the XSE- 
LECT package and the average spectrum was extracted from the summed event file. Events 
for the spectral analysis were selected within a circle with a 20 pixel (~ 47arcsec) radius, 
which encloses about 95% of the PSF, centered on the source position. The background was 
extracted from a nearby circular region of 40 pixel radius. The source spectrum was binned 
to ensure a minimum of 20 counts per bin to utilize the \ 2 minimization fitting technique. 
In addition, we needed to apply a small energy offset (~ 40 eV) to the observed 
energy spectrum. The origin of this correction is likely to be CCD charge traps 
generated by radiation and high-energy proton damage (SWIFT-XRT-CALDB- 
12), which affect mostly the lowest energies (first 1-2 bins) of the spectrum. 
The ancillary response files were generated with the xrtmkarf task applying corrections for 
the PSF losses and CCD defects using the cumulative exposure map. The latest response 
matrices (v.011) available in the Swift CALDB were used. 

The XRT average spectrum in the 0.3 — 10 keV energy band was fitted using the XSPEC 
package. We adopted a log-parabolic model of the form F{E) = K ■ (^y)~ ( - r+l3 ' l ° 9 ^'^^ 



( iMassaro et al.l l2004al lbl) with an absorption hydrogen-equivalent column d ensity fixed to 
the G alactic value in the direction of the source, which is 1.61 x 10 20 cm -2 (IKalberla et al. 



20051 ). We found that this model provided a good description of the observed spectrum, with 



the exception of the 1.4—2.3 keV energy band where spectral fit residuals were present. These 



residuals are due to known XRT calibration uncertainties (SWIFT-XRT-CALDB-12 125 !) and 
hence we decided to exclude the 1.4 — 2.3 keV energy band from the analysis. The resulting 
spectral fit gave the following parameters: K = (1.839 ± 0.002) x 10 _1 ph cm -2 s _1 keV -1 , 
T = 2.178 ± 0.002, = 0.391 ± 0.004. The XRT SED data shown in figure [8] were corrected 
for the Galactic absorption and then binned in 16 energy intervals. 



5.1.5. RXTE/PCA 



The Rossi X-ray Timing Explorer (RXTE; iBradt et al.l 119931 ) satellite performed 59 
pointing observations of Mrk421 during the time interval MJD 54851 and 54972. These 



http : / /heasarc . gsf c . nasa.gov/docs/heasarc/ caldb/ swift /docs/xrt/SWIFT-XRT-CALDB-09_v!2 .pdf 
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observations amount to a total exposure of 118 ks, which was requested through a dedicated 
Cycle 13 proposal to provide X-ray coverage for this multi-instrument campaign on Mrk421. 

The data analysis was performed using FTOOLS v6.9 and following the procedures and 



filtering criteria recommended by the RXTE Guest Observer Facility 126 ! after September 
2007. The average net count rate from Mrk421 was about 25 ct s -1 per pcu (in 
the energy range 3 — 20 keV) with flux variations typically much smaller than 
a factor of two. Consequently, the observations were filtered following the conservative 
procedures for faint sources: Earth elevation angle greater than 10°, pointing offset less than 
0.02°, time since the peak of the last SAA (South Atlantic Anomaly) passage greater than 
30 minutes, and electron contamination less than 0.1. For further details on the analysis 



of faint sources with RXTE, see the online Cook Booh 127 !. In the data analysis, in order to 
increase the quality of the signal, only the first xenon layer of PCU2 was used. We used the 
package pcabackest to model the background and the package saextrct to produce spectra 



for the source and background files and the script! I pcarsp to produce the response matrix. 



The PCA average spectrum in the 3 — 32 keV energy band was fitted using the XSPEC 
package using a PL function with an exponential cutoff (cutoffpl) with a non- variable neutral 
Hydrogen column density Nh fixed to the Galactic value in the direction of the source 



(1.61 x 10 20 cm~ 2 ; (IKalberla et al.l 120051 )). However, since the PCA bandpass starts at 3 
keV, the value for Nh used does not significantly affect our results. The resulting spectral 
fit provided a good representation of the data for the following parameters: normalization 
parameter K = (2.77±0.03) x 10" 1 ph cm" 2 s" 1 keV~\ photon index T = 2.413 ±0.015, and 
cutoff energy E exp = 22.9 ± 1.3 keV. The obtained 23 energy bins PCA average spectrum is 
shown in Figure 



5.1.6. Swift/ BAT 



The Swift Burst Alert Telescope (BAT: lBarthelmy et al.ll2005l ) analysis results presented 
in this paper were derived with all the available data during the time inter val MJD 548 5 and 



54983. The spectrum was extracted following the recipes presented in lAjello et al.l (12008 



2009_bJ). This spectrum is constructed by weight averaging the source spectra extracted over 
short exposures (e.g. 300 s) and it is representative of the averaged source emission over the 



12( 



http : //www . universe . nasa . gov/xrays/ programs/rxte/pca/doc/bkg/bkg-2007-saa/ 



12 



http : / /heasarc . gsf c . nasa . gov/docs/xte/recipes/ cook_book .html 



128 



The CALDB files are located at http://heasarc.gsfc.nasa.gov/FTP/caldb 
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time range spanned by t he observations. Th ese spectra are accurate to the mCrab level and 
the reader is referred to lAjello et al.l ( )2009al ) for more details. The Swift /BAT spectrum in 
the 15-200 keV energy range is consistent with a PL function with normalization parameter 
K = 0.46 ± 0.27 ph cm~ 2 s" 1 keV" 1 and photon index T = 3.0 ± 0.3. The last two 
flux points are within one standard deviation from the above mentioned PL 
function and hence the apparent upturn given by these last two data points in 
the spectrum is not significant. 



5.1.7. MAGIC 



MAGIC is a system of two 17m-diameter IACTs for VHE 7-ray astronomy located on 
the Canary Island of La Palma, at an altitude of 2200 m above sea level. At the time of the 
observation, MAGIC-II, the new second telescope of the current array system, was still in 
its commissioning phase so that Mrk421 w as observed in stan d-alone mode by MAGIC-I, 
which is in scientific operation since 2004 ([Albert et al.l 120081) . The MAGIC observations 
were performed in the so-called "wobble" mode ( lDaumlll997l ). In order to have a low energy 
threshold, only observations at zenith angles less than 35° were used in this analysis. The bad 
weather and a shut down for a scheduled hardware system upgrade during the period MJD 
54948-54960 (April 27-May 13) significantly reduced the amount of time that had initially 
be en scheduled f or thi s ca mpaign. The data were analyzed following the prescription given 
by lAlbert et al.l (120081 ) and lAliu et al.l ( 120091 ). The data surviving the quality cuts amounted 
to a total of 27.7 hours. The preliminary reconstructed photon fluxes for the individual 
observations gave an average flux of about 50% that of the Crab Nebula, with relatively mild 
(typically less than factor 2) flux variations. The derived spectrum was unfold ed to correct 
for the effects of the limited energy resolution of the detector and possible bias ( 1 Albert et al. 
2007a ). The resulting spectrum was fit satisfactorily with a single log-parabola function: 
F(E) = K ■ ( 3 ^ eV )~( r+ ^ og (o.3Tcv)). The resulting spectral fit gave the following parameters: 
K = (6.50 ± 0.13) x 10~ 10 ph cm" 2 s" 1 erg" 1 , T = 2.48 ± 0.03, (3 = 0.33 ± 0.06, with 
X 2 /NDF = 11/6. A fit with a simple power-law function gives ^/NDF = 47/7, which 
confirmed the existence of curvature in the VHE spectrum. 



5.2. Fermi-hAT Spectra During the Campaign 

The Mrk421 spectrum measured by Fermi-LAT during the period covered by the mul- 
tifrequency campaign is shown in panel (b) of Fig[71 The spectrum can be described with 
a single PL function with photon index 1.75 ± 0.03 and photon flux F(> 0.3 GeV) = 
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(6.1 ± 0.3) x 10 8 phcm 2 s l ; which is somewhat lower than the average spectrum over the 
first 1.5 years of Fermi-LAT operation (see Figure[5]). 

For comparison purposes we also computed the spectra for the time periods before and 
after the multifrequency campaign (the time intervals MJD 54683-54850 and MJD 54983- 
55248, respectively). These two spectra are shown in panels (a) and (c) of Figured The 
two spectra can be described very satisfactorily with single PL functions of photon indices 
1.79 ±0.03 and 1.78 ±0.02 and photon fluxes F(> 0.3 GeV) = (7.1 ±0.3) x 10~ 8 phones -1 
and F(> 0.3 GeV) = (7.9 ± 0.2) x 10~ 8 ph cm _2 s± Therefore, during the multifrequency 
campaign, Mrk421 showed a spectral shape that is compatible with the periods before and 
after the campaign, and a photon flux which is about 20% lower than before the campaign 
and 30% lower than after the campaign. 



5.3. The Average Broad Band SED during the Multifrequency Campaign 

The average SED of Mrk421 resulting from our 4.5-month-long multifrequency cam- 
paign is shown in Figure[HJ This is the most complete SED ever collected for Mrk421 or 
for any oth e r BL L ac object (although an SED of nearly similar quality was reported in 



Abdo et al.l (l2011bl ) for Mrk 501). At the highest energies, the combination of Fermi-LAT 
and MAGIC allows us to measure, for the first time, the high energy bump without any 
gap; both the rising and falling segments of the components are precisely determined by 
the data. The low energy bump is also very well measured; Swift /BAT and RXTE/PCA 
describe its falling part, Swift/XRT the peak, and the Swift/UV and the various optical 
and IR observations describe the rising part. The rising tail of this peak was also measured 
with various radio instruments. Especially important are the observations from SMA at 
225 GHz, which help connecting the the bottom (radio) to the peak (optical/X-rays) of the 
synchrotron bump (in the vF v representation). The flux measurements by VLBA, especially 
the ones corresponding to the core, provide us with the radio flux density from a region 
that is presumably not much larger than the blazar emission region. Therefore, the radio 
flux densities from interferometric observations (from the VLBA core) are expected to be 
close upper limits to the radio continuum of the blazar emission component. On the other 
hand, the low frequency radio observations performed with single dish instruments have a 
relatively large contamination from the non-blazar emission and are probably considerably 
above the energy flux from the blazar emission region. The only spectral intervals lack- 
ing observations are from 1 meV - 0.4eV, and 200 keV - 100 MeV, where the 
sensitivity of the current instruments is insufficient to detect Mrk 421. We note 
however, that the detailed GeV coverage together with our broadband, 1-zone 
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SSC modeling strongly constrains the expected emission in the difficult to access 
1 meV - 0.4 eV bandpass. 

During this campaign, Mrk421 showed low activity and relatively small flux 
variations at all frequencies (|Panequell2009l ) . At VHE (> 100 GeV), the measured 
flux is half the flux from the Crab N ebula, which is am o ng the lowest fluxe s 
recorded by MAGIC for this source ( lAlbert et all l2007al : lAleksic et all bold ). 
At X-rays, the fluxes observed during this campaign are about 15 mCrab, which 
is about 3 times higher than the lowest fluxes measured by RXTE/ASM since 
1996. Therefore, because of the low flux, low (multifrequency) variability and 
the large density of observations, the collected data during this campaign can 
be considered an excellent proxy for the low/quiescent state SED of Mrk421. 
It is worth stressing that the good agreement in the overlapping energies of the various 
instruments (which had somewhat different time coverage during the campaign) supports 
this hypothesis. 



6. SED Modeling 



We turn now to modeling the multifrequency data set collected during the 
4.5 month campaign in the context of homogeneous hadronic and leptonic mod- 
els. The models discussed below assume emission mainly from a single, spheri- 
cal and homogeneous region of the jet. This is a good approximation to model 
flaring events with observed correlated variability (where the dynamical time 
scale does not exceed the flaring time scale significantly), although it is an over- 
simplification for quiescent states, where the measured blazar emission might be 
produced by the radiation from different zones characterized by different values 
of the relevant parameters. There are several models in the literature along 
those lines (e.g. iGhisellini et al.l l2005l: iKatarzvhski et al.ll2008k iGraff et al1 l2008: 



Giannios et al 



2009) but at the cost of introducing more free parameters that 



are, consequently, less well constrained and more difficult to compare between 
models. This is particularly problematic if a "limited" data set (in a time and 
energy coverage) is employed in the modeling, although it could work well if the 
amount of multifrequency data is extensive enough to substantially constrain 
the parameter space. In this work we adopted the 1-zone homogeneous models 
for their simplicity as well as for being able to compare with previous works. 
The 1-zone homogeneous models are the most widely used models to describe 
the SED of high-peaked BL Lacs. Furthermore, although the modeled SED is 
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averaged over 4.5 months of observations, the very low observed multifrequency 
variability during this campaign, and in particular the lack of strong keV and 
GeV variability (see Figs [1] and ED in these timescales suggests that the presented 
data are a good representation of the average broad-band emission of Mrk421 
on timescales of few days. We, therefore, feel confident that the physical param- 
eters required by our modeling to reproduce the average 4.5 month SED are a 
good representation of the physical conditions at the emission region down to 
timescales of a few days, which is comparable to the dynamical timescale derived 
from the models we discuss. The implications (and caveats) of the modeling re- 
sults are discussed in section 

Mrk421 is at a relatively low redshift (z = 0.031), yet the attenuation of its VHE 
MAGIC spectrum by the extragalactic background light (EBL) is non-negligible for all mod- 
els and hence needs to be accounted for using a parameterization for the EBL density. The 
EBL absorption at 4 TeV, the highest energy bin of the MAGIC data (absorption will be 
less at lower energ ies), varies accordin g to the model used from e~ Tl1 = 0.29 for the "Fast 
Evolu tion" model (ptecker et al. 120061 ) to e~ T " n = 0.58 for the models of iFranceschini et al. 



(12008( 1 and 



Gilmore et al 



model of Finkeetal 



(120091 ). with most models giving e 



0.5-0.6, in cluding the 



J2OI0J) and the "best fit" model of iKneiske et all (l2004h. We have 



de-absorbed the TeV data from MAGIC with the IFranceschini et al.l (120081 ) model, although 
most other models give comparable results. 



6.1. Hadronic Model 

If relativistic protons are present in the jet of Mrk 421, hadronic interactions, if above the 
interaction threshold, must be considered for modeling the source emissi on. For the present 



mode ling we use the hadronic Synchrotron-Proton Blazar (SPB) model of lMiicke et al.l ( 12001 



20031 ). Here, the relativistic electrons (e) injected in the strongly magnetized (with homoge- 
neous magnetic field with strength B) blob lose energy predominantly through synchrotron 
emission. The resulting synchrotron radiation of the primary e component dominates the 
low energy bump of the blazar SED, and serves as target photon field for interactions with 
the instantaneously injected relativistic protons (with index a p = a e ) and pair (synchrotron- 
supported) cascading. 

Figures [9] and [TO] show a satisfactory (single zone) SPB model representation of the 
data from Mrk 421 collected during the campaign. The corresponding parameter values 
are reported in Table [3j In order to fit the optical data, the lowest energy of the injected 
electrons is required to be maintained as 7 e min ~ 700 through the steady state. This requires 
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a continuous electron injection rate density of at least > 1.4cm _3 s 1 to balance 
the synchrotron losses at that energy, and is about a factor ~ 100 larger than the proton 
injection rate. The radio fluxes predicted by the model are significantly below the observed 
8-230 GHz radio fluxes. This is related to the model being designed to follow the evolution 
of the jet emission during 7-ray production where radiative cooling dominates over adiabatic 
cooling. Here, the emission region is optically thick up to ~100 GHz frequencies, and the 
synchrotron cooling break (j e ~ 10) would be below the synchrotron-self- absorption turn- 
over. The introduction of additional, poorly constrained components would be necessary 
to account for the subsequent evolution of the jet through the expansion phase where the 
synchrotron radiation becomes gradually optically thin at cm wavelengths. This is omitted 
in the following modeling. 

The measured spectra in the 7-ray band (> 1 GeV) is dominated by synchrotron radi- 
ation from short-lived muons (produced during photomeson production) as well as proton 
synchrotron radiation, with significant overall reprocessing, while below this energy the tt- 
cascade dominates. The interplay between muon and proton synchrotron radiation together 
with appreciable cascade synchrotron radiation initiated by the pairs and high energy pho- 
tons from photomeson production, is responsible for the observed MeV-GeV flux. The TeV 
emission is dominated by the high energy photons from the muon synchrotron component. 
The source intrinsic model SED predicts > 10 TeV emission on a level of 2 to 3 orders of 
magnitude below the sub- TeV flux, which, will be further weakened by 7-ray absorption by 
the EBL. 

The overall required particle and field energy density are within a factor 5 
of equipartition, and a total jet power (as measured in the galaxy rest frame) of 
4 x 10 44 erg s 1 in a greement with expectations for a weakly accreting disk of a 



BL Lac object (see ICaoll2003h . 



Alternative model fits are possible if the injected electron and proton components do 
not have the same power- law index. This "relaxation" of the model would add one extra 
parameter and so would allow for improvement in the data-model agreement, especially 
around the synchrotron peak and the high energy bump. It would also allow a larger tolerance 
on the size region R, which is considered to be small in the SPB model fit presented here. 



6.2. Leptonic Model 



The simplest leptonic model typically used to describe the emission from BL Lac ob- 
jects is the 1-zone Synchrotron Self-Compton model (SSC). Within this framework, the radio 
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through X-ray emission is produced by synchrotron radiation from electrons in a homoge- 
neous, randomly-oriented magnetic field (B) and the the 7-rays are produced by inverse 
Compton scattering of the synchrotron photons by the sa me electrons which produce them. 
For this purpose, we use the 1-zone SSC code described in lFinke et al.l ( 120081 ). The electron 
distribution from 1-zone SSC models is typically parameterized with one or two power-law 
(PL) functions (that is, zero or one break) within an electron Lorentz factor range defined by 
7 min and ■jmax (where the electron energy is •ym e c 2 ). We use the same approach in this work. 
However, we find that, in order to properly describe the shape of the measured broad-band 
SED during the 4.5 months long campaign, the model requires an electron distribution pa- 
rameterized with three PL functions (and hence two breaks). In other words, we must add 2 
extra free parameters to the model: the second break at 7^,2 and the index of the third PL 
function p 3 . Note that a second break was also needed to d escribe the SED of Mrk 501 in the 
context of the synchrotron/SSC model (lAbdo et al.ll2011bl ). An alternative possibility might 
be to use an electron distribution para metrized with a curve d function such as that result- 
ing fro m episodic partic l e acce leration ( iPerlman et al.l 120051 ) or the log-parabolic function 
used in lTramacere et al.l (120091 ) . However, we note that such a parameterization might have 
problems describing the highest X-ray energies, where the current SED data (RXTE/PCA 
and Swift /BAT) do not show a large spectrum curvature. 

Even though the very complete SED constrains the shape of the electron distribution 
quite well, there is still some degeneracy in the range of allowed values for the general source 
parameters R (comoving blob radius), B and 5 (doppler factor). For a given break in the 
measured low energy (synchrotron) bump, the break in the electron distribution 7^ scales 
as 1/ V B5. In order to minimize the range of possible parameters, we note that the emitting 
region radius is constrained by the variability time, t v , so that 

Set ftct 



ii, = <~ 

1 + z ~ 1 + z 



1) 



During the observing campaign, Mrk 421 was in a rather low activi ty state, with mul- 
tifrequency flux variations occurring on timescales larger than 1 day ( jPanequd 120091 ) . so 
we used t vm i n = 1 day in our modeling. In addition, given that this only gives an up- 
per limit on the size scale, and the hi story of fast variability detected for this object (e.g. 
Gaidos et al.lll996l ; iGiebels et al.l 120071 ). we also performed the SED model using t Vtmin = 1 
hour. The resulting SED models obtained with these two variability timescales are shown in 
Figure [TTJ with the parameter values reported in table HI The blob radii are large enough in 
these models that synchrotron self-absorption (SSA) is not important; for the t VjTn i n = 1 hr 
model, vssa = 3 x 10 10 Hz, at which frequency a break is barely visible in Figure [TlJ It is 
worth stressing the good agreement between the model and the data: the model describes 
very satisfactorily the entire measured broad-band SED. The model goes through the SMA 
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(225 GHz) data point, as well as through the VLBA (43 GHz) data point for the partially- 
resolved radio core. The size of the VLBA core of the 2009 data from Mrk421 at 15 GHz 
and 43 GHz is ~ 0.06-0.12 mas (as reported in Section I5.1.ip or, using the conversion scale 
0.61pc/mas ~ 1-2 xl0 17 cm. The VLBA size estimation is the FWHM of a Gaussian repre- 
senting the brightness distribution of the blob, which coul d be approximated as 0.9 times the 
radius of a corresponding spherical blob (IMarscherl 1 19831 ) . That implies that the size of the 
VLBA core is comparable (factor ~2-4 larger) than that of the model blob for t var = 1 day 
(~ 5 x 10 16 cm). Therefore, it is reasonable to consider that the radio flux density from the 
VLBA core is indeed dominated by the radio flux density of the blazar emission. The other 
radio observations are single dish measurements and hence integrate over a region that is 
orders of magnitude larger than the blazar emission. Consequently, we treat them as upper 
limits for the model. 

The powers of the different jet components derived from the model fits (as- 
suming r = 8) are also reported in Table HJ, Estimates for the mass o f the super 



mass i ve black hole in Mrk421 range from 2 x 10 8 M Q to 9 x 10 s M (IBarth et al. 



20031 ; IWu et al.ll2002l ). and hence the Eddington luminosity should be between 
2.6 x 10 46 and 1.2 x 10 47 erg s _1 , that is well above the jet luminosity. 

It is important to note that the parameters resulting from the modeling of our broad- 
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)09|). One difference, as already noted, is that an extra break is required. This could 
be a feature of Mrk 421 in all states, but we only now have the simultaneous high quality 
spectral coverage to identify it. For the model with t var = 1 day (which is the time variability 
observed during the multifrequency campaign), additional differences with previous models 
are in R, which is an order of magnitude larger, and B, which is an order of magnitude 
smaller. This mostly results from the longer variability time in this low state. Note that 
using a shorter variability (t var = 1 hour, green curve) gives a smaller R and bigger B than 
most models of this source. 

Another difference in our 1-zone SSC model with respect to previous works relate to 
the parameter ■jmin- This parameter has typically not been well constrained because the 
single-dish radio data can only be used as upper limits for the radio flux from the blazar 
emission. This means that the obtained value for 7 m j n (for a given set of other parameters 
R, B, and 5) can only be taken as a lower limit: a higher value of 7 m i„ is usually possible. In 
our modeling we use simultaneous Fermi-LAT data as well as SMA and VLBA 
radio data, which we assume is dominated by the blazar emission. We note that 
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the size of the emission from our SED model fit (when using t var ~1 day) is 
comparable to the partially resolved VLBA radio core and hence we think this 
assumption is reasonable. The requirement that the model SED fit goes through 
those radio points puts further constrains to the model, and in particular to the 
parameter 7 mm : a decrease in the value of 7 mm would over-predict the radio data, 
while an increase of 7 m j„ would under-predict the SMA and VLBA core radio 
data, as well as the Fermi-LAT spectrum below 1 GeV if the increase in 7 min 
would be large. We explored model fits with different 7 mm and pi, and found 
that, for the SSC model fit with t var = 1 day (red curve in Figure [TT]) . j m i n is well 
constrained within a factor of 2 to the value of 8 x 10 2 (see Figure 1121) . In the 
case of the SSC model with t var = 1 hour (green curve in Figure [TT]) . if we make 
the same assumption that the SMA and VLBA core emission is dominated by 



the blazer emission 129 !. 7 min can be between 2 x 10 2 up to 10 3 , and still provide a 
good match to the SMA/ VLB A/optical data and the Fermi-LAT spectrum. In 
any case, for any variability timescale, the electron distribution does not extend 
down to 7 m i n ~ 1 to a few, and is constrained within a factor of 2. This is 
particularly relevant because, for power-law distributions with index p > 2, the 
jet power carried by the electrons is dominated by the low energy electrons. 
Therefore, the tight constraints on 7 mm translate into tight constraints on the jet power 
carried by the electrons. For instance, in the case of the model with t var = 1 hour, using 
7mm = 10 3 (instead of 7 mm = 4 x 10 2 ) would reduce the jet power carried by electrons from 
Pj >e « 10 44 erg s' 1 down to Pj i€ « 8 x 10 43 erg s^ 1 . 

Another parameter where the results presented here differ from previous results in the 
literature is the first PL index p\. This parameter is dominated by the optical and UV data 
points connecting with the Swift /XRT, as well as by the necessity of matching the model 
with the Fermi-LAT GeV data. Note that our model fit also goes over the SMA and VLBA 
(partially-resolved) core fluxes. Again, since these constrains did not exist (or were not used) 
in the past, most of the 1-zone SSC model results (for Mrk421) in the literature report a pi 
value that differs from the one reported in this work. We note however that the values for the 
parameters pi and p^ from our model fits, which are constrained mostly by the X-ray/TeV 
data, are actually quite similar to the parameters p\ and p 2 from the previous 1-zone SSC 
model fits to Mrk421 data. 



129 In the case of t var ~ 1 hour, the size of the emission region derived from the SSC model is one order of 
magnitude smaller than the size of the VLBA core and hence the used assumption is somewhat less valid 
than for the model with t var ~ 1 day 



-30 - 



7. Discussion 

In this section of the paper we discuss the implications of the experimental 
and SED modeling results presented in the previous sections. As explained at 
the beginning of section 6, for simplicity and for the sake of comparison with 
previously published results, we modeled the SED with scenarios based on 1- 
zone homogeneous emitting regions, which are commonly used to parameterize 
the broad-band emission in blazars. We note that this is a simplification of the 
problem; the emission in blazar jets could be produced in an inhomogeneous or 
stratified region, as well as in N independent regions. An alternative and quite 
realistic scenario could be a standing shock where particle acceleration takes 
place a nd radiation is being produced as the jet flow or superluminal knots cross 



it (e.g. iKomissarov &: Falld Il997l ; iMarscher et al.ll2008l ). The Lorentz factor of 



the plasma, as it flows through the standing (and by necessity oblique) shock 
is the Lorentz factor (and through setting the angle, the Doppler factor) of the 



model. We note however that, as discussed in ISikora et al.l (119971 ). the steady- 
sate emission could also be parametrized by N moving blobs that only radiate 
when passing through the standing shock. If at any given moment, only one of 
these blobs were visible at the observer frame, the 1-zone homogeneous model 
could be a plausible approximation of the standing-shock scenario. 

In any case, the important thing is that, in the proposed physical scenario, 
the stability time-scale of the particle accelerating shock front is not connected 
to the much shorter cooling times that give rise to spectral features. For as long 
as the injection of particles in the blob and the dynamics of the blob remain 
unchanged, the SED, along with the breaks due to radiative cooling and due 
to the value of 7 m j n where Fermi acceleration presumably picks up, will remain 
unchanged. The lack of (substantial) multifrequency variability observed during 
this campaign suggests that this is the case, and hence that the 4.5-months- 
averaged SED is also representative of the broad-band emission of SED during 
much shorter periods of time that are comparable to the dynamical timescales 
derived from the models. 



7.1. What are the Spectral Breaks Telling Us? 

In our homogenous leptonic model we reproduce the location of the vj v peaks by fitting 
the Lorentz factors and ( a s we U a s the values of B and 8) where the electron 

energy distribution breaks. There is, however, a Lorentz factor where one typically (in 
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blazar modeling) expects a break in the electron energy distribution (EED), and this is the 
Lorentz factor 7 C = 37rm e c 2 / (a T B 2 R) where the escape time from the source equals the 
radiative (synchrotron) cooling time. The fact that the values of the second break, jbrk,2, 
fit by our leptonic models (7^2 = 3.9 x 10 5 , 1.7 x 10 5 ) are similar to the Lorentz factors 
(7 C = 1.6 x 10 5 , 3.3 x 10 5 ) where a cooling break in the EED is expected, strongly suggests 
that the second break in the EED derived from the modeling is indeed the cooling break. 

The observed spectral shape in both the low and high energy SED components are 
reproduced in our homogenous model by a change of electron index Ap = p 3 — p 2 = 2.0. 
Such a large break in the EED is in contrast to the canonical cooling break Ap = p-$—p2 = 1.0 
that produces a spectral index change of Aa = 0.5, as predicted for homogenous models (e.g., 
Longair 1994). An attempt to model the data fixing Ap = P3—P2 = 1-0 gave unsatisfactory 
results, and hence this is not an option; a large spectral break is needed. It would be tempting 
to speculate that what we observe is not a cooling break, but rather something that results 
from a characteristic of the acceleration process which is not understood and that, therefore, 
does not bind us to the Ap = 1.0 constraint. But we would then have to attribute to shear 
fortuity the fact that the Lorentz factors where this break takes place are very close to the 
Lorentz factors where cooling is actually expected. 

The question that naturally arises is why, although the EED break postulated by the 
homogeneous model is at nearly the same energy as the expected cooling break, the spectral 
break observed is stronger. Such strong breaks are the rule rather than t he exception in some 
non-thermal sources like pulsar-wind nebulae and extragalactic jets (see iReynoldsl 120091 ) and 
the explanations that have been given relax the assumption of a homo genous emitting zone, 
invoking gradients in the physical quantities describing the system ( iMarscherl Il980[ ). In 
all inhomogenous models, electrons are injected at an inlet and are advected downstream 
suffering radiative losses that result in the effective size of the source declining with increasing 
frequency for a given spectral component. In sources where the beaming of the emitted 
radiation is the same throughout the source (this is the case for non-relativistic flows or for 
relativistic flows with small velocity gradients), the spectral break formed is stronger than 
the canonical Aa = 0.5 if the physical conditions change in such a way that the emissivity at 
a given frequency increases downstream (Wilson 1975, Coleman & Bicknell 1988, Reynolds 
2009). 



If, in addition to these considerations, we a 



dients, either in the form of a dece lerating flow (IGeorganopoulos et al.l 120031 ) or the form of 



low for significant relativist ic velocity gra 



a fast spine and slow sheath flow ( iGhisellini et al.l 120051 ). the resulting differential beaming 
of the emitted radiation can result in spectral breaks stronger that Aa « 0.5. Studies of 
the SEDs of sources with different jet orientations (e.g. radio galaxies and blazars) can help 
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to understand the importance of differential beaming, and therefore of relativistic velocity 
gradients in these flows. Because in all these models the volume of the source emitting at a 
given frequency is connected to the predicted spectral break, it should be possible to use the 
variability time-scale at different frequencies to constrain the physics of the inhomogenous 
flow. 



7.2. Physical Properties of Mrk 421 

As mentioned in §5.31 the SED emerging from the multifrequency campaign is the most 
complete and accurate representation of the low/quiescent state of Mrk 421 to date. This 
data provided us with an unprecedented opportunity to constrain and tune state-of-the-art 
modeling codes. In §|6]we modeled the SED within two different frameworks: a leptonic and 
a hadronic scenario. Both models are able to represent the overall SED. As can be seen in 
Figures M and [LTJ the leptonic model fits describe the observational data somewhat better 
than the hadronic model; yet we also note here that, in this paper, the leptonic model has one 
more free parameter than the hadronic model. A very efficient way of discriminating between 
the two scenarios would be through multi-wavelength variability observations. It is however 
interesting to discuss the differences between the two model descriptions we presented above. 



7.2.1. Size and location of the emitting region 



The characteristic size to which the size of the emitting region must be compared is the 
gravitational radius of the Mrk 42 1 black hole. For a black hole mass of ~ 2 — 9 x 10 8 M Q 



( IBarth et al.ll2003t IWu et al.ll2002[ ). the corresponding size is R g « 0.5 — 2.0 x 10 cm. In 



the hadronic model the source size can be as small as R = 4 x 10 14 cm (larger source sizes 
can not be ruled out though; see Sect. 6.1), within one order of magnitude of the gravita- 
tional radius. The consequence is a dense synchrotron photon energy density that facilitates 
frequent interactions with relativistic protons, resulting in a strong reprocessed/cascade com- 
ponent which leads to a softening of the spectrum occurring mostly below 100 MeV. The 
Fermi-LAT analysis presented in this paper (which used the instrument response function 
given by P6_V3_DIFFUSE) is not sensitive to these low energies and hence the evaluation of 
this potential softening in the spectrum will have to be done with future analyses (and more 
data). This will potentially allow accurate determination of spectra down to photon energies 
of ~ 20 MeV with the LAT. 



The leptonic model can accommodate a large range of values for R, as long as it is not so 
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compact that internal 77 attenuation becomes too strong and absorbs the TeV 7-rays. In the 
particular case of t var = 1 day, which is supported by the low activity and low multifrequency 
variability observed during the campaign, R = 5 x 10 16 cm, that is 2-3 orders of magnitude 
larger than the gravitational radius. Under the assumption that the emission comes from 
the entire (or large fraction of the) cross-section of the jet, and assuming a conical jet, the 
location of the emitting region would be given by L ~ R/@, where 9 ~ 1/T ~ 1/8. Therefore, 
under these assumptions, which are valid for large distances (L 3> R g ) when the outflow is 
fully formed, the leptonic model would put the emission region at L ~ 10 3 — 10 A R g . We note 
however that, since the R for the leptonic model is considered an upper limit on the blob 
size scale (see eqn. [T]), this distance should also be considered as an upper limit as well. 



7.2.2. Particle content and particle acceleration 

The particle content predicted by the hadronic and leptonic scenarios are different by 
construction. In the hadronic scenario presented in §6.1[ the dominantly radiating parti- 
cles are protons, secondary electron/positron pairs, muons, and pions, in addition to the 
primary electrons. In the leptonic scenario, the dominantly radiating particles are the pri- 
mary electrons only. In both cases, the distribution of particles are clearly non-thermal and 
acceleration mechanisms are required. 

In the leptonic scenario, the PL index p± = 2.2, which is the canonical particle spectral 
index from efficient lst-order Fermi acceleration at the fronts of relativistic shocks, suggests 
that this process is at work in Mrk421. For electrons to be picked up by lst-order Fermi 
acceleration in perpendicular shocks, their Larmor radius is required to be significantly larger 
than the width of the shock, which for electron-proton plasmas is set by the Larmor radius 
of the dynamically dominant particles (electrons or protons). The large 7 m j„ (= 8 x 10 2 ) 
provided by the model implies that electrons are efficiently accelerated by the 
Fermi mechanism only above this energy and that below this energy they are 
accelerated by a different mechanism that produces an extremely hard electron 
dis tribution. Suc h pre-acceleration mechanisms have been discussed in the past (e.g 



• 1 



Hoshino et al.lll992l ). The suggestion that the Fermi mechanism picks up only after 7 min 
(= 8 x 10 2 ) suggests a large thickness of the shock, which would imply that the shock is 
dominated by (cold ) protons. We refer the reader to the Fermi-LAT paper on Mrk501 



(lAbdo et al.ll2011bl ) for more detailed discussion on this topic. In addition, in §6.21 and §7.11 
we argued that the second break r ybrk,2fn e c 2 (~ 200 GeV) is probably due to synchrotron 
cooling (the electrons radiate most of their energy before existing the region of size R), but 
the first break •jbrk,im e c 2 (~ 25 GeV ) must be related to the acceleration mechanism; and 
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hence the leptonic model also requires that electrons above the first break are accelerated less 
efficiently. At this point it is interesting to note that the 1-zone SSC model of Mrk501 in 2009 
(where the source was also observed mostly in a quiescent state), returned ihrk ^m P c 2 ~ 2 



GeV w ith essentially the same spectral change (0.5) in the electron distribution (lAbdo et al. 



201 lbl ). Therefore, the first break (presumedly related to the acceleration mechanism) is of 
the same magnitude and located approximately at the same energy for both Mrk421 and 
Mrk501, which might suggest a common property in the quiescent state of HSP BL Lac 
objects detected at TeV energies. 

The presence of intrinsic high energy breaks in the EE D electron energy distrib ution 



has been o bserved in several o f the Fermi-L AT blazars (see lAbdo et al.l l2009al |2010aj) . As 



reported in lAbdo et al.l feOlOal ). this characteristic was observed on several FSRQs, and it is 
present in some LSP-BLLacs, and a small n umber of ISP-BLLa cs; yet it is absent in all 1LAC 
HSP-BLLacs. In this paper (as well as in lAbdo et al.l (j2011bl )) we claim that such feature 
is also present in HSP-BLLacs like Mrk421 and Mrk501, yet for those objects, the breaks 
in the EED can only be accessed through proper SED modeling because they are smaller 
in magnitude, and somewhat smoothed in the high energy component. We note that, for 
HSP-BLLacs, the high energy bump is believed to be produced by the EED upscattering 
seed-photons from a wide energy range (the synchrotron photons emitted by the EED itself) 
and hence all the features from the EED are smoothed out. On the other hand, in the 
other blazar objects like FSRQs, the high energy bump is believed to be produced by the 
EED upscattering (external) seed-photons which have a "relatively narrow" energy range. 
In this later case (external compton), the features of the EED may be directly seen in 
the gamma-ray spectrum. Another interesting observation is that, at least for one of the 
FSRQs, 3C 454 .3, the location and th e magnitude of the break seems to be insensitive to 
flux variations (lAckermann et al.l 120101 ) . If the break observed in Mrk421 and Mrk501 is of 
the same nature than that of 3C 454.3, we should also expect to see this break at the same 
location (~ 20 GeV) regardless of the activity level of these sources. 

In the hadronic scenario of Figure [91 the blazar emission comes from a compact (R ~ 
a few R g ) highly magnetized emission region, which should be sufficiently far away from 
the central engine so that the photon density from the accretion disk is much smaller than 
the density of synchrotron photons. The gyroradius of the highest energy protons (Rl = 
1 P ,maxmpC 2 / (eB) in Gaussian-cgs units) is ~ 1.4 x 10 14 cm, which is a factor of ~3 times 
smaller than the radius of the spherical region responsible for the blazar emission (R = 
4 x 10 14 cm), hence (barely) fulfilling the Hillas criterium. The small size of the emitting 
region, the ultra-high particle energies and the somewhat higher (by factor ~5) particle 
energy density with respect to the magnetic energy density imply that this scenario requires 
extreme acceleration and confinement conditions. 
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7.2.3. Energetics of the jet 

The power of the various components of the flow differ in the two models. In the SPB 
model, the particle energy density is about a factor of ~5 higher than the magnetic field 
energy density and the proton energy density dominates over that of the electrons by a 
factor of ~ 40. In the leptonic model the electron energy density dominates over that of 
the magnetic field by a factor of 10. By construction, the leptonic model does not constrain 
the proton content and hence we need to make assumptions on the number of protons. It is 
reasonable to use charge neutrality to justify a comparable number of electrons and protons. 
Under this assumption, the leptonic model predicts that the energy carried by the electrons 
(which is dominated by the parameter 7 m j n ~ 10 3 ) is comparable to that carried by the 
(cold) protons. 

The overall jet power determined by the hadronic model is Pj e t = 4.4 x 10 44 erg s _1 . 
For the day variability timescale leptonic model, assuming one cold proton per radiating 
electron, the power carried by the protons would be 4.4 x 10 43 erg s _1 , giving a total jet 
power of Pj e t = 1.9 x 10 erg s . In both cases, the computed jet power is a small fraction 
(~ 10~ 2 — 10~ 3 ) of the Eddington luminosity for the supermassive black hole in Mrk421 
(2 • 1O 8 M ) which is L Edd ~ 10 46 - 10 47 erg s" 1 . 



7.3. Interpretation of the Reported Variability 

In §0 we reported the 7-ray flux/spectral variations of Mrk421 as measured by the 
Fermi-LAT instrument during the first 1.5 years of operation. The flux and spectral index 
were determined on 7-day-long time intervals. We showed that, while the 7-ray flux above 
0.3 GeV flux changed by a factor of ~ 3, the PL photon index variations are consistent with 
statistical fluctuations (Figure [Q and the spectral variability could only be detected when 
comparing the variability in the 7-ray flux above 2 GeV with that one from the 7-ray flux 
below 2 GeV. It is worth pointing out that, in the case of the TeV blazar Mrk501, the 7-ray 
flux above 2 GeV was also found to vary more than the 7-ray flux below 2 GeV. Yet unlike 
Mrk421, Mrk501 was less bright at 7-rays and the flux variations above 2 GeV seem to be 
larger, which produced statistically si gnificant changes in the photon index from the PL fit 



in the energy range 0.3-400 GeV (see lAbdo et al.ll2011bl ). In any case, it is interesting to 



note that in these two (classical) TeV objects, the flux variations above few GeV are larger 
than the ones below a few GeV, which might suggest that this is a common property in HSP 
BL Lac objects detected at TeV energies. 



In §3] we also showed (see Figures El [31 and H]) that the X-ray variability is significantly 
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higher than that in the 7-ray band measured by Fermi-LAT . In addition, we also saw that 
the 15 — 50keV (BAT) and the 2 — lOkeV (ASM) fluxes are positively correlated, and that 
the BAT flux is more variable than the ASM flux. In other words, when the source flares in 
X-rays, the X-ray spectrum becomes harder. 

In order to understand this long baseline X-ray/7-ray variability within our leptonic 
scenario, we decomposed the 7-ray bump of the SED into the various contributions from the 
various segments of the EED, according to our 1-zone SSC model, in a similar way as it was 



done in iTavecchio et al.l ( 119981 ) . This is depicted in Figure [TSJ The contributions of different 
segments of the EED are indicated by different colors. As shown, the low-energy electrons, 
Imin < 7 < 7&r,i, which are emitting synchrotron photons up to the observed frequencies 
of ~ 5.2 x 10 15 Hz, dominate the production of 7-rays up to the observed photon energies 
of ~20 GeV (green line). The contribution of higher energy electrons with Lorentz factors 
76r,i < 7 < Ibr, 2 is pronounced within the observed synchrotron range 5 x 10 15 — 10 17 Hz, 
and at 7-ray energies from ~20 GeV up to ~TeV (blue line). Finally, the highest energy 
tail of the electron energy distribution, 7 > 7^2, responsible for the production of the 
observed X-ray synchrotron continuum (> 0.5 keV) generates the bulk of 7-rays with the 
observed energies > TeV (purple line). Because of the electrons upscattering the broad energy 
range of synchrotron photons, the emission of the different electron segments are somewhat 
connected, as shown in the bottom plot of Figure [13j Specifically, the low energy electrons 
have also contributed to the TeV photon flux through the emitted synchrotron photons which 
are being upscattered by the high energy electrons. Hence, changes in the number of low 
energy electrons should also have an impact on the TeV photon flux. However, note that 
the synchrotron photons emitted by the high energy electrons, which are up-scattered in the 
Klein-Nishina regime, do not have any significant contribution to the gamma-ray flux, thus 
changes in the number of high energy electrons (say 7 > 7^2) will not significantly change 
the MeV/GeV photon flux. 

Within our 1-zone SSC scenario, the 7-rays measured by Fermi-LAT are mostly pro- 
duced by the low energy electrons (7 < 7b r , 1) while the X-rays seen by ASM and BAT are 
mostly produced by the highest energy electrons (7 > 7^,2)- In this scenario, the signifi- 
cantly higher variability in the X-rays with respect to that of 7-rays suggests that the flux 
variations in Mrk421 are dominated by changes in the number of the highest energy elec- 
trons. Note that the same trend is observed in the X-rays (ASM versus BAT) and 7-rays 
(below versus above 2 GeV); the variability in the emission increases with the energy of the 
radiating electrons. 



The greater variability in the radiation produced by the highest energy electrons is not 
surprising. The cooling timescales of the electrons from synchrotron and inverse Compton 
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(in the Thomson regime) losses scale as t oc 7 _1 , and hence it is expected that the emission 
from higher energy electrons will be the most variable. However, since the high energy 
electrons are the ones losing their energy fastest, in order to keep the source emitting in X- 
rays, injection (acceleration) of electrons up to the highest energies is needed. This injection 
(acceleration) of high energy electrons could well be the origin of the flux variations in 
Mrk421. The details of this high energy electron injection could be parameterized by changes 
in the parameters 7^2, P3 and 7 ma x within the framework of the 1-zone SSC model that 



could result from episodic acceleration events ( jPerlman et al.ll2005l ). The characterization of 
the SED evolution (and hence SSC parameter variations) will be one of the prime subjects of 
the forthcoming publications with the multi-instrument variability and correlation during the 
campaign s in 2009 130 l and 2 01Q 131 I. We note here that SSC models, both one-zone and multi- 



zone (e.g. iGraff et al.ll2008l ). predict a positive correlation between the X-rays and the TeV 
7-rays measured by IACTs. Indeed, during the 2010 campaign the source was detected in a 
flaring state with the TeV instruments (see ATel #2443). Such an X-r ay/TeV correlation has 
been established in the past for this object (see iMaraschi et al.lll999l ). although the relatio n 
is not simple. Sometimes it is linear and other times quadratic (e.g. iFossati et al.ll2008l ). 
The complexity of this correlation is also consistent with our 1-zone SSC model; the X-rays 
are produced by electrons with 7 > 7^2, while the TeV photons are produced by electrons 
with 7 > 7br,i, and is indirectly affected by the electrons with 7 < 7^1 through the emitted 
synchrotron photons that are used as seed photons for the inverse Compton scattering (see 
the bottom plot of Figure [T3|) . 



We also note that the 1-zone SSC scenario presented here predicts a direct correlation on 
the basis of simultaneous data sets between the low energy gamma-rays (from Fermi) and the 
sub-millimeter (SMA) and optical frequencies, since both energy bands are produced by the 
lowest energy electrons in the source. On the other hand, our SPB model fit does not require 
such a strict correlation, but there could be a loose correlation if electrons and protons are 
accelerated together. In particular, a direct correlation with zero time lag between the mm 
radio frequencies and the 7-rays is not expected in our SPB model because the radiation 
at these two energy bands are produced at different sites. The radiation in the X- and 
7-ray band originates from the primary electrons, and from the protons and secondary 
particles created by proton-initiated processes, respectively. Consequently, although a loose 
correlation between the X-ray and 7-ray band can be expected if protons and electrons are 



130 For details of the 2009 campaign, see the URL 

https : //confluence . slac . Stanford . edu/display/GLAMC0G/Campaign+on+Mrk421+ ( Jan+2009+to+May+2009) 

131 p or details of the 2010 campaign, see the URL 

https : //confluence . slac . Stanford. edu/display/GLAMCQG/Campaign+on+Mrk421+'/,28December+2009+to+December+2010'/,29 
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accelerated together, a strict correlation with zero time lag is rather unlikely in our model 
fit. 

During the 2009 and 2010 campaigns, Mrk421 was very densely sampled during a very 
long baseline (4.5 and 6 months for the 2009 and 2010 campaign, respectively) and hence 
these data sets will provide excellent information for performing a very detailed study of 
these multi-band relations. In particular, during the campaign in 2010, there were regular 
observations with VLBA and SMA, which will allow us to study with a greater level of detail 
the relationship between the rising parts of the low energy and high energy bumps, where 
the predictions from the leptonic and hadronic models differ. 

8. Conclusions 

In this work, we reported on the 7-ray activity of Mrk421 as measured by the LAT 
instrument on board the Fermi satellite during its first 1.5 years of operation, from 2008 
August 5 (MJD 54683) to 2009 March 12 (MJD 55248). Because of the large leap in capa- 
bilities of LAT with respect to its predecessor, EGRET, this is the most extensive study of 
the 7-ray activity of this object at GeV photon energies to date. The Fermi-LAT spectrum 
(quantified with a single power-law function) was evaluated for 7-day-long time intervals. 
The average photon flux above 0.3 GeV was found to be (7.23 ± 0.16) x 10~ 8 ph cm -2 s _1 , 
and the average photon index 1.78 ± 0.02. The photon flux changed significantly (up to a 
factor ~ 3) while the spectral variations were mild. The variations in the PL photon index 
were not statistically significant, yet the light curves and variability quantification below and 
above 2 GeV showed that the high 7-ray energies vary more than the low energy 7-rays. We 
found F var (E < 2 GeV) = 0.16 ± 0.04 while F var (E > 2 GeV) = 0.33 ± 0.04. We compared 
the LAT 7-ray activity in these two energy ranges (0.2 — 2 GeV and > 2 GeV) with the 
X-ray activity recorded by the all-sky instruments RXTE/ASM (2- lOkeV) and Swift/BAT 
(15 — 50keV). We found that X-rays are significantly more variable than 7-rays, with no 
significant (<2<r) correlation between them. We also found that, within the X-ray and 7-ray 
energy bands, the variability increased with photon energy. The physical interpretation of 
this result within the context of the 1-zone SSC model is that the variability in the radiation 
increases with the energy of the electrons that produce them, which is expected given the 
radiating time scales for synchrotron and inverse Compton emission. 

We also presented the first results from the 4. 5- month- long multifrequency campaign 
on Mrk421, which lasted from 2009 January 19 (MJD 54850) to 2009 June 1 (MJD 54983). 
During this time period, the source was systematically observed from radio to TeV energies. 
Because of the low activity and low variability shown during this campaign, the compiled 
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data provided us with the best SED yet of Mrk421 in the low/quiescent state. 

The broadband SED was modeled with two different scenarios: a leptonic (1 zone SSC) 
and a hadronic model (SPB). Both frameworks are able to describe reasonably well the 
average SED, implying comparable powers for the jet emission, which constitute only a 
small fraction (~ 10~ 2 — 10 -3 ) of the Eddington luminosity. However, those models differ on 
the predicted environment for the blazar emission: the leptonic scenario constrains the size 
to be R < 10 4 R g , the magnetic field to B ~ 0.05 G and particles (electrons) with energies 
up to ~ 5 • 10 13 eV while, if a e = a p , our hadronic scenario implies a size of the emitting 
region of a few R g , a magnetic field B ~ 50 G and particles (protons) with energies up to 
~ 2 • 10 18 eV, which requires extreme conditions for particle acceleration and confinement. 

The leptonic scenario suggests that the acceleration of the radiating particles (electrons) 
is through diffusive shock acceleration in relativistic shocks mediated by cold protons, and 
that this mechanism accelerates particles (electron s) less efficiently abo ve an energy of ~ 25 



GeV, which is comparable to what was reported in lAbdo et al.l ( 1201 lbl ) for another classical 
TeV blazar, Mrk501. In addition, unlike what was observed for Mrk501, in the case of 
Mrk421 a stronger-than-canonical electron cooling break was required to reproduce the 
observed SED, which might suggest that the blazar emitting region is inhomogeneous. 

Within the SSC model (Figure ITT]) , the observed X-ray/7-ray variability during the 
first 1.5 years of Fermi operation indicates that the flux variations in Mrk421 are produced 
by acceleration of the highest energy electrons, which radiate in the X-ray and TeV bands, 
and lose energy, radiating as they do so in the optical and GeV range. In our hadronic 
model (Figure [9]), a rather loose correlation between the X- and 7-ray band is expected 
if electrons and protons are accelerated together. A forthcoming publication will report on 
whether these emission models can reproduce the multi-band flux variations observed during 
the intensive campaigns on Mrk421 performed in 2009 and 2010. Those studies should help 
us distinguish between the hadronic and the leptonic scenarios and eventually lead to a 
better understanding of one of the fundamental mysteries of blazars: how flux variations are 
produced. 
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Fig. 1. — Left: 7-ray flux at photon energies above 0.3 GeV (top) and spectral photon index 
from a power-law fit (bottom) for Mrk421 for 7-day-long time intervals from 2008 August 
5 (MJD 54683) to 2009 March 12 (MJD 55248). Vertical bars denote la uncertainties and 
the horizontal error bars denote the width of the time interval. The black dashed line and 
legend show the results from a constant fit to the entire data set. Right: scatter plot of the 
photon index versus flux. 



Table 1. Discrete Correlation Function (DCF) computed using the flux values reported in 
Figure [2] (7-day-long time intervals, first 1.5 years of Fermi operation) and Figure [3] 
(3-day- long time interval during the last 5 months, where the X-ray activity was high). 

The DCF values are given for time lag zero. 



Interval 


ASM-BAT 


ASM-LAT<2 GeV 


ASM-i AT> 2 GeV 


BAT-LAT<2 GeV 


BAT-LAT>2 GeV 


LAT K 2GeV — LAT >2 GeV 


7-day-long 


0.73 ± 0.20 


0.28 ±0.15 


0.35 ± 0.14 


0.20 ±0.13 


0.26 ±0.13 


0.31 ±0.14 


3-day-long 


0.65 ± 0.13 


0.01 ±0.18 


0.15 ± 0.19 


-0.03 ± 0.13 


0.01 ±0.13 


0.29 ±0.17 



Note. — The DCF was computed as prescribed in lEdelson fc Krolilj dl98Sh . 
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Fig. 2. — Multifrequency light curves of Mrk421 with 7-day-long time bins obtained with 
3 all-sky-monitoring instruments: RXTE/ASM (2 - lOkeV, first); Swift /BAT (15 - 50keV, 
second) and Fermi-LAT for two different energy ranges (0.2 — 2 GeV, third, and > 2GeV, 
fourth). The light curves cover the period from 2008 August 5 (MJD 54683) to 2009 March 
12 (MJD 55248). Vertical bars denote la uncertainties and horizontal error bars show the 
width of the time interval. The black dashed lines and legends show the results from constant 
fits to the entire data set. The vertical dashed lines denote the time intervals with 
the extensive multifrequency campaigns during the 2009 and 2010 seasons. 
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Fig. 3. — Multifrequency light curves of Mrk421 with 3-day-long time bins obtained with 3 
all-sky-monitoring instruments: RXTE/ASM (2 - lOkeV, first); Sw/i/ BAT (15 - 50keV, 
second) and Fermi-LAT for two different energy ranges (0.2 — 2 GeV, third, and > 2GeV, 
fourth). The light curves cover the period from 2009 October 4 to 2010 March 12. Vertical 
bars denote la uncertainties and horizontal error bars show the width of the time interval. 
The black dashed lines and legends show the results from constant fits to the entire data set. 
The vertical dashed lines denote the beginning of the extensive multifrequency 
campaign on Mrk421 during the 2010 season. 
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Fig. 4. — Fractional variability parameter for 1.5 year data (2008 August 5- 2009 March 
12) from 3 all-sky-monitoring instruments: RXTE/ASM (2-10 keV, first); Swift/BAT (15- 
50 keV, second) and Fermi-LAT for 2 en ergy ranges 0.2-2 Ge V and 2-300 GeV. The fractional 
variability was computed according to IVaughan et al.l ( 120031 ) using the light curves from 
Figure [2j Vertical bars denote la uncertainties and horizontal bars indicate the width of 
each energy bin. 
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Fig. 5. — Fermi spectrum of Mrk421 during the period from 2008 August 5 to 2010 February 
20. Black line is the likelihood PL fit; red contour is the 68% uncertainty of the fit and the 
black data points show the energy fluxes computed on differential energy ranges. The inlay 
summarizes the unbinned likelihood PL fit in the energy range 0.3-400 GeV. 
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Fig. 6. — Time and Energy coverage during the multifrequency campaign. For the sake of 
clarity, the minimum observing time displayed in the plot was set to half a day. 



Table 2. List of instruments participating in the multifrequency campaign and used in the 

compillation of the SED shown in FigureE] 



Instrument/Observatory 



Energy range covered 



Web page 



MAGIC 

Whipple" 

Fermi-LAT 

Swift /BAT 

RXTE/PCA 

Swift /XRT 

Swift/VYOT 

Abastumani (through GASP-WEBT program) 
Lulin (through GASP-WEBT program) 

Roque de los Muchachos (KVA) (through GASP-WEBT program) 

St. Petersburg (through GASP-WEBT program) 

Talmassons (through GASP-WEBT program) 

Valle d'Aosta (through GASP-WEBT program) 

GRT 

ROVOR 

New Mexico Skies 

MitSume 

OAGH 

WIRO 

SMA 

VLBA 

Noto 

Metsahovi (through GASP-WEBT program) 
VLBA (through MOJAVE program) 
OVRO 
Mcdicina 

UMRAO (through GASP-WEBT program) 
RATAN-600 

Effelsberg (through F-GAMMA program) 



0.08-5.0 TeV 
0.4-2.0 TeV 
0.1-400 GeV 
14-195 keV 
3-32 keV 
0.3-9.6 keV 

UVW1, UVM2, UVW2 

R band 

R band 

R band 

R band 

R band 

R band 

V, R, B, I bands 
B, R, V bands 
R, V bands 
g, Rc, Ic bands 
H, J, K bands 
J, K bands 
225 GHz 

4.8, 8.3, 15.4, 23.8, 43.2 GHz 

8.4, 22.3 GHz 

37 GHz 

15 GHz 

15 GHz 

8.4 GHz 

4.8, 8.0, 14.5 GHz 

2.3, 4.8, 7.7, 11.1, 22.2 GHz 

2.6, 4.6, 7.8, 10.3, 13.6, 21.7, 31 GHz 



http : //wwwmagic . mppmu . mpg . de/ 



http : //veritas . sao . arizona. edu/content/blogsection/6/40/ 
: //www-glast . Stanford. edu/index .html 
: //heasarc .gsfc .nasa.gov/docs/swif t/swiftsc .html/ 
: //heasarc . gsfc .nasa. gov/docs/xte/rxte .html 
: //heasarc .gsfc .nasa. gov/docs/swif t/swiftsc .html 



http 



: //heasarc .gsfc .nasa. gov/docs/swif t/swiftsc .html 
: //www. oato . inaf . it/blazars/webt/ 
: //www. oato . inaf . it/blazars/webt/ 
: //www. oato . inaf . it/blazars/webt/ 
: //www. oato . inaf . it/blazars/webt/ 
: //www. oato . inaf . it/blazars/webt/ 
: //www. oato . inaf . it/blazars/webt/ 



K1 



http: //asd. gsfc .nasa.gov/Takanori . Sakamoto/GRT/index r html 

^ 1 ' 

tp : //rovor . byu . edu/ 

tp: //www.nmskies . com/equipment .html/ 



http : //www . vlba . nrao . edu/ 



: //www . hp . phys . t itech . ac . jp/mitsume/index . html 
: //astro . inaoep .mx/en/observatories/oagh/ 
: //physics . uwyo . edu/- chip/wiro/wiro . html 
: //smal . sma.hawaii . edu/ 



http://www.noto . ira. inaf . it/ 

http://www.metsahovi.fi/ 

http : //www . physics . purdue . edu/MQ JAVE/ 

http: //www. ovro . caltech. edu/ 

http: //www. med. ira. inaf . it/index_EN .htm 

http: //www. oato . inaf . it/blazars/webt/ 



http: //w0 . sao . ru/ratan/ 



http: //www. mpif r-bonn. mpg. de/div/ef f elsberg/index_e .html/ 



Note. — The energy range shown in column 2 is the actual energy range covered during the Mrk421 observations, and not the instrument nominal energy range, which might 
only be achievable for bright sources and excellent observing conditions. 



Note. - 



(a) The Whipple spectra were not included in Figure[8] See text for further comments. 
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Fig. 7. — Fermi spectra of Mrk421 for several time intervals of interest. The panel (a) shows 
the spectrum for the time period before the multifrequency campaign (MJD 54683-54850), 
the panel (b) for the time interval corresponding to the multifrequency campaign (MJD 
54850-54983) and the panel (c) for the period after the campaign (MJD54983-55248). In 
all panels, the black line depicts the result of the unbinned likelihood PL fit and the red 
contours denote the 68% uncertainty of the PL fit. The legend reports the results from the 
unbinned likelihood PL fit in the energy range 0.3 — 400 GeV. 
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Fig. 8. — Spectral energy distribution of Mrk421 averaged over all the observations taken 
during the multifrequency campaign from 2009 January 19 (MJD 54850) to 2009 June 1 
(MJD 54983). The legend reports the correspondence between the instruments and the mea- 
sured fluxes. The host galaxy has been subtracted, and the optical/X-ray data were corrected 
for the Galactic extinction. The TeV data from MAGIC were corr ected for the absorption i n 
the extragalactic background light using the prescription given in iFranceschini et al.l ( 120081 ). 
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Fig. 9.— Hadronic model fit components: 7r°-cascade (black dotted line), n^ 1 cascade 
(green dashed-dotted line), /i-synchrotron and cascade (blue dashed-triple-dotted line), pro- 
ton synchrotron and cascade (red dashed line). The black thick solid line is the sum of all 
emission components (which also includes the synchrotron emission of the primary electrons 
at optical/X-ray frequencies). The resulting model parameters are reported in table El 
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Fig. 10. — Expanded view of the high energy bump of the SED data and model presented 
in Figure M 
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Table 3. Parameter values from the SPB model fit to the SED from Mrk421 shown in 

Figure EJ 



Parameter 


Symbol 


Value 


Doppler Factor 


S 


12 


Magnetic Field [G] 


B 


50 


Comoving blob radius [cm] 


R 


4 x 10 14 


Power law Index of the injected electron distribution a 


a e 


1.9 


Power law Index of the injected proton distribution a 


a p 


1.9 


Minimum Electron Lorentz Factor 


'Je,min 


7 x 10 2 


Maximum Electron Lorentz Factor 


ye, max 


4 x 10 4 


Minimum Proton Lorentz Factor b 


'Jp^min 


1 


Maximum Proton Lorentz Factor 


yp^max 


2.3 x 10 9 


Energy density in protons [erg cm -3 ] 


Ul p 


510 


Ratio of number of electrons with respect to protons 


e/p 


90 


Jet Power [erg s ] 


Pjet 


4.5 x 10 44 



Note. — (a) The model assumes a e = a p , hence only one free parameter. 

Note. — (6) The parameter 7 Pjm i n was fixed to the lowest possible value, 1, and 
hence this is actually not a free parameter. 
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Fig. 11. — SED of Mrk421 with two 1-zone SSC model fits obtained with different minimum 
variability timescales: t var = 1 day (red curve) and t var = 1 hour (green curve) . The 
parameter values are reported in Table HJ See text for further details. 
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Table 4. Parameter values from the 1-zone SSC model fits to the SED from Mrk421 

shown in Figure [TTJ 



Parameter 


Symbol 


Red 


Curve 


Green Curve 


Variability Timescale [s] a 




8.64 x 10 4 


3.6 x 10 3 


Doppler Factor 


5 




21 


50 


Magnetic Field [G] 


B 


3.8 


x 10~ 2 


8.2 x 10~ 2 


Comoving blob radius [cm] 


R 


5.2 


x 10 16 


5.3 x 10 15 


Low-Energy Electron Spectral Index 


Pi 




2.2 


2.2 


Medium-Energy Electron Spectral Index 


P2 




2.7 


2.7 


High-Energy Electron Spectral Index 


P3 




4.7 


4.7 


Minimum Electron Lorentz Factor 


Tmin 


8.0 


x 10 2 


4 x 10 2 


Breakl Electron Lorentz Factor 


Ibrkl 


5.0 


x 10 4 


2.2 x 10 4 


Break2 Electron Lorentz Factor 


7brk2 


3.9 


x 10 5 


1.7 x 10 5 


Maximum Electron Lorentz Factor 


'Jmax 


1.0 


x 10 8 


1.0 x 10 8 


Jet Power in Magnetic Field [erg s~ 1 ] b 


P 3 ,B 


1.3 


x 10 43 


3.6 x 10 42 


Jet Power in Electrons [erg s — 1 ] 


P*,e 


1.3 


x 10 44 


1.0 x 10 44 


Jet Power in Photons [erg s -1 ] 6 




6.3 


x 10 42 


1.1 x 10 42 



Note. — (a) The variability timescale was not derived from the model fit, but rather 
used as an input (constrain) to the model. See text for further details. 

Note. — (b) The quantities Pj t B and Pj :V h are derived quantities; only Pj >e is a free 
parameter in the model. 



-60 - 




Fig. 12. — SSC model fit of the SED from Mrk421 presented in Figure [IT] (for t var ~1 
day), with variations by a factor of 2 of the parameter 7™„, together with adjustments in 
the parameter pi in order to match the experimental data. See text for further details. 
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Fig. 13. — Decomposition of the high energy bump of the SSC continuum for Mrk421. 
The data points are the same as in the high energy bump from Figure [TTJ The SSC fit 
to the average spectrum is denoted by the red solid curve. Top: Contributions of the 
different segments of electrons comptonizing the whole synchrotron continuum (green curve: 
7mm < 7 < lbr,u blue curve: 7^1 < 7 < 7^2; purple curve: 7^2 < 7)- Bottom: 
Contributions of the different segments of electrons (as in the top panel) Comptonizing 
different segments of the synchrotron continuum (solid curves: v < Vbr,i — 5.3 x 10 15 Hz; 
dashed curves: z/ fcr 1 < v < v hr 2 — 1.3 x 10 17 Hz; dotted curves, corresponding to v > ^,2) 



